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Abstract 

In this paper we establish lower and upper Gaussian bounds for the solutions to the 
heat and wave equations driven by an additive Gaussian noise, using the techniques of 
Malliavin calculus and recent density estimates obtained by Nourdin and Viens in [19J. 
In particular, we deal with the one-dimensional stochastic heat equation in [0,1] driven 
by the space-time white noise, and the stochastic heat and wave equations in M'' (d > 1 
and d < 3, respectively) driven by a Gaussian noise which is white in time and has a 
general spatially homogeneous correlation. 



1 Introduction 

The Malliavin calculus, or stochastic calculus of variations, is a suitable technique for proving 
that a given random vector F = (Fi,...,Fm) on the Wiener space possesses a smooth 
probability density. There has been a current interest in the applications of the stochastic 
calculus of variations to obtain lower and upper bounds of Gaussian type for the density of 
a given Wiener functional. The starting point of this research is the work by Kusuoka and 
Stroock [12J, where they proved that the density of a uniformly hypoelliptic diffusion whose 
drift is a smooth combination of its diffusion coefficient has a lower bound of Gaussian type. 
Recently, three different approaches have been developed to derive Gaussian-type bounds for 
densities of general Wiener functional using the techniques of Malliavin calculus: 
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(i) In the paper [TT], Kohatsu-Higa proposes a general methodology for computing lower 
bounds for a multidimensional functional F = (Fi, . . . , Fm) of a Wiener sheet in [0, T] x 
W^. The method uses an approximation of F by means of a sequence of conditionally 
non-degenerate random variables adapted to the filtration generated by the white noise. 
This paper was inspired by the work by Kusuoka and Stroock [T2], and it can be used 
in a non Markovian framework. As an application, the author obtains lower bounds 
for the solution to the one-dimensional heat equation driven by a space-time white 
noise, assuming that the diffusion coefficient a is bounded away from zero. The ideas 
introduced in this paper were later developed in the work by Bally pLj to get Gaussian 
lower bounds for locally elliptic Ito processes. 

(ii) Nourdin and Viens in [19J have proved a new formula for the density of a one-dimensional 
Wiener functional in terms of the Malliavin calculus (see (12. 2 ft below). As an application, 
they obtain upper and lower Gaussian bounds for the density of the maximum of a general 
Gaussian process. 

(iii) In [14j Malliavin and Nualart derived Gaussian lower bounds for multidimensional Wiener 
functionals under an exponential moment condition on the divergence of a covering 
vector field. A one-dimensional version of this result was obtained by Nualart in [TU] . 

The purpose of this paper is to apply the results obtained by Nourdin and Viens in [IS] 
to the solutions of different classes of stochastic partial differential equations driven by an 
additive Gaussian noise. Upper and lower Gaussian estimates for the density of solutions to 
stochastic partial differential equations are essential tools in the potential analysis for these 
type of equations (see the recent works [6, 7, 8J). 

The paper is organized as follows. In Section 2 we recall briefly the results of [T^]. Section 
3 deals with the one-dimensional heat equation on the interval [0, 1] with Dirichlet boundary 
conditions and driven by a space-time white noise. Let us denote by u{t,x) its solution 
evaluated at some {t,x) E IR+ x (0,1) and set m = E\u{t,x)\. Then, we derive a gaussian 
lower and upper bound for the density of u{t, x) of the form 



The lower bound differs from that obtained by Kohatsu-Hida in the paper [TT], where the first 
factor was Ct* . 

Section 3 is devoted to the stochastic heat equation in with an additive Gaussian noise 
which is white in time and it has a spatially homogeneous correlation in space. We assume 
that the spectral measure of the noise integrates (1 + |(^P)'' for some r] E (0, 1), and we obtain 
lower and upper Gaussian bounds for the density of the solution of the form (11.11) . that involve 
the powers t and t^^"^ (see Theorem 14.41 for the precise statement). Let us mention at this 
point that it has been of much importance the fact the Malliavin derivative of the solution 
defines a non-negative process, since it solves a linear parabolic equation. 

Finally, in Section 4 we consider a stochastic wave equation in dimension d = 1,2,3, again 
driven by an additive Gaussian noise which is white in time and it has a spatially homogeneous 
correlation in space. In this case we obtain lower and upper bounds only for t small enough, and 
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they involve the powers and t^~'^'^, where again rj is given by the integrability of the spectral 
measure. We have not been able to overcome this technical restriction on the time parameter 
because, in comparison with the heat equation, the Malliavin derivative of the solution in this 
case does not need to be a non-negative function neither a function itself. 

2 Gaussian density estimates 

in this section, we recall a general method set up in in order to show that a smooth 
random variable in the sense of Malliavin calculus has a probability density admitting Gaussian 
lower and upper bounds. 

First of all, let us briefly describe the Gaussian context in which we will be working on 
and the main elements of the Malliavin calculus associated to it. Namely, suppose that 
in a complete probability space P) we are given a centered Gaussian family W = 

{W{h), h E H} of random variables indexed by a separable Hilbert space H with covariance 

E{W{h)W{g)) = {h,g)H, KgeH. 

The family W is usually called an isonormal Gaussian process on H. Assume that JF is the 
a-field generated by W. 

Let us use standard notation for the main operators of the Malliavin calculus determined 
by the family W (see, for instance, |20]). More precisely, we denote by D the Malliavin 
derivative, defined as a closed and unbounded operator from L^(r2) into L'^{Vt]H), whose 
domain is denoted by D^'^. The adjoint of the operator D is denoted by 5, usually called the 
divergence operator. A random element u E L'^iVt; H) belongs to the domain of 5 if and only 
if it satisfies 

\E{{DF,5{u))H)\<C\\F\\L2(^n). 

for any F E D^'^, where the constant C only depends on u. For any element u in the domain 
of 6, the random variable 5{u) can be characterized by the duality relationship 

E{F6{u)) = E{{DF,u)h), 

for every F E D^-^. 

Any random variable F in L'^{n, J^, P) can be decomposed by means of its Wiener chaos 
expansion (see [2U1 Section 1.1.1]), which is usually written as 

oo 

F = Y,JnF, 

n=0 

where J„ denotes the projection onto the nth Wiener chaos. 

Using the chaos expansion one may define the operator L by the formula L = Yl'^=o ~^'^n, 
which is called the generator of the Orstein-Uhlenbeck semigroup (see [201 Section 1.4]). It is 
related to the Malliavin derivative D and its adjoint S through the formula 

6DF = -LF, 



3 



in the sense that F belongs to the domain of L if and only if it belongs to the domain of 5D, 
and in this case the above equality holds. 

One can also define the inverse of L, denoted by L^^, as follows: for any F E L^(i7, JF, P), 
set L-^F := T.n=i-^JnF. Then it holds that LL'^F = F - E{F), for any F e 
L^(r2,jF, P), so that L^^ acts as the inverse of L for centered random variables. 

Let us consider F E D^'^ with mean zero and we define the following function in R: 

g{z) := E[{DF, -DL-'F)h\F = z]. 

By [18, Proposition 3.9], it holds that g(z) > on the support of F. Then, [19, Theorem 
3.1] states that if the random variable g{F) is bounded away from zero almost surely, that is 

9{F) > ci > 0, a.s., 

for some constant Ci, then F has a density p whose support is M satisfying, almost everywhere: 

As a consequence, see JW, Corollary 3.4], if one also has that g{F) < C2, a.s., then the density 
p satisfies, for almost all 2; G M: 

^l^lexpf-ilUpW,i|^e.pf-ilV (2.3) 



2ci V 2c2y ~ 2c2 V 2c_ 

Let us also mention that [121 Proposition 3.5] provides a more suitable formula for g{F) 
for computational purposes. Indeed, given a random variable F E D^'^, one can write DF = 
where is a measurable mapping from to H, determined P o iy~^-almost 
surely (see l2Qj, p. 54-55). Then 

POO 

g{F)= / e-^E {^F{W),^F{e-^W + \/i^^^W'))h\F dd, (2.4) 
Jo L -I 

where W stands for an independent copy of W such that W and W are defined on the 
product probability space {Vt x n',J^0j^',P x P'). Eventually, E denotes the mathematical 
expectation with respect to P x P'. 

Formula (12.41) can be still rewritten in the following form: 



g{F) = J^ e-^E E' {JdF.DF^ ^\F de, (2.5) 



where, for any random variable X defined in {Vt, ^, P), X denotes the shifted random variable 

in ^7 X f2' 

X{iu,u') = X(e-^c<;+ v^r^e^w')> uj E u' E n. 

Notice that, indeed, X depends on the parameter 9, but we have decided to drop its explicit 
dependence for the sake of simplicity. 

Along the paper we denote by C a generic constant which may vary from line to line. 
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3 Density estimates for the stochastic heat equation in 

[0,1] 

In this section, we will consider a stochastic heat equation in [0,1], with Dirichlet boundary 
conditions, some non-linear drift b and with an additive space-time white noise perturbation. 
We aim to give sufficient conditions on the coefficient b ensuring Gaussian lower and upper 
bounds for the probability density of the solution at any point. 

3.1 The stochastic heat equation in [0,1] 

We are concerned with the following one dimensional heat equation driven by a space-time 
white noise: 

— (t, x) - —{t, x) = b{u{t, x)) + aW{t, x), {t, x) e [0, T] x [0, 1], (3.6) 

where T > 0, the initial condition is given by a continuous function uq : [0, 1] ^ M and we 
consider Dirichlet boundary conditions. That is, 

u{0,x) = uo{x), X G [0,1], (3.7) 

u{t,0) = u{t, 1) = 0, te [0,T]. 

The real-valued random field solution to Equation (13. 6 ft will be denoted by {u{t,x), {t,x) E 
[0,T] X [0,1]}. The function 6 : M — M is of class having a bounded derivative and 
cr > is a constant. We assume that {W{t,x), {t,x) G [0,T] x [0, 1]} is a Brownian sheet 
on [0,T] X [0,1], defined in a complete probability space (i7,jF, P). That is, {W{t,x)} is a 
centered Gaussian family with the covariance function 

E{W{t, x)W{s, y)) = {t A s){x Ay). 

For < t < T, let JF^ be the cr-field generated by the random variables {W{s 
[0,t] X [0,1]} and the P-null sets. 

The solution to the formal Equation (13.61 ) is understood in the mild sense: a {jFt}-adapted 
stochastic process {u{t,x), {t,x) G [0,T] x [0,1]} solves (13.61 ) with initial and boundary 
conditions (EJ]) if, for any {t,x) G (0,T] x (0,1), 

u{t,x)= Gt{x,y)uo{y)dy + / Gt-s{x,y)b{u{s,y))dyds 
Jo Jo Jo 

' Gt^,{x,y)W{ds,dy), (3.8) 



a 




JO 



where Gt{x,y), {t,x,y) G IR+ x (0,1)^, denotes the Green function associated to the heat 
equation on [0, 1] with Dirichlet boundary conditions. We will use the following facts, for any 

< t' < t and X G (0,1), 

0<Gt{x,y)<^e-^, (3.9) 
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c,Vt^<f f \Gt-s{x,y)\^dyds<^Vt^, (3.10) 
Jt' Jo vStt 

where is a positive constant. We also have mfa<x<i-aCx > for any a E (0,1). In the 
case of Neumann boundary conditions, Cx does not depend on x (see, for instance, [T71 Lemma 
A1.2]). 

Let us also mention that the stochastic integral in (13.80 is understood as an integral with 
respect to the Brownian sheet in the sense of Walsh [26J. 

Existence and uniqueness of mild solution for Equation (13.81) can be deduced from the 
results in [26]; for an even more general setting, see also [2j. Malliavin calculus applied to 
Equation (13.81 ) have been dealt with in [3J. In this case, we consider the Gaussian context 
associated to the space-time white noise. That is, we have H = L^(]R+ x [0,1]) and the 
Gaussian family {W{h), h E L'^{R+ x [0, 1])} is given by the Wiener integral 

W{h)= [ [ h{s,y)Wids,dy). 
Jr+ Jo 

A consequence of Proposition 4.3 and Theorem 2.2 in [Tj is that, for all {t,x) E (0,T] x 
[0,1], the random variable u(t,x) belongs to D^'^, the Malliavin derivative satisfies the linear 
parabolic equation 

Dr,zu{t, x) = aGt^rix, ^) + J Gt-s{x, y)h'{u{s, y))Dr,zu{s, y) dyds, (3.11) 

for any (r, z) E [0, t] x [0, 1] and {t, x) E [0, T] x [0, 1] , and the probability law of u{t, x) has 
a density. The positivity of a and G guarantees that the solution of Equation (13. lift remains 
non-negative, that is Dr,zu{t,x) > 0, a.s. This will be a key point in proving that the density 
of u{t,x) admits Gaussian upper and lower estimates. 



3.2 Gaussian bounds for the density 

We fix T > and consider {u{t,x), {t,x) E [0,T] x [0,1]} the unique mild solution of 
Equation (13.8ft . This section will be devoted to proof the following result. 

Theorem 3.1 Assume that the drift coefficient b is a function with bounded derivative. 
Then, for all t E (0,T] and x E (0,1), the random variable u{t,x) possesses a density p 
satisfying the following statement: for almost every z eM., 

E\u(t,x)—m\ f (z — mY) , , E\u{t,x^—m\ [ (z — m)?^ 

' ^' ^ 'exp ^ / \<v{z)< ' ^' \ ^exp ^ i . (3-12) 



C2U I Clt2 J Cit2 [ Git 

where m := E{u{t, x)) and Gi, G2 are positive quantities depending on a, ||6'||oo, T and x. 

The statement of the above Theorem 13.21 will be a consequence of formula (12.21) (see 
[T^ Theorem 3.1]) and the following proposition. In order not to overload notations, set 
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F := u{t,x) — E{u{t,x)), and recall that, as it has been specified in Section |2] (see (12.5 
therein), 



gin 



e-'E 



E' 




'0 ^0 



DrzF [DrzF] dzdr 



F 



de 



F' yj J F)r,zu{t,x) (jJr,zu{t,x)^ dzdr 



de, 



(3.13) 



where DF = {DF){e'^uj + VT^^e^u'). 

Proposition 3.2 Fix T > and assume that the drift coefficient b is of class and has a 
bounded derivative. There exist positive constants Ci, C2 such that 

Cit^ < 9{F) < C2t^ , (3.14) 

for all te (o,r]. 

In order to prove Proposition 13. 2[ we will need the following technical lemma: 

Lemma 3.3 Let t > 0. Assume that b E and it has a bounded derivative. There exists a 
positive constant K depending on a and Hf^'Hoo- such that, for any 6 E (0, 1]; 



sup 

(i-s)t<v<t J n-S)t Jo 
0<?y<l 



1 

E[\Dr,M^,y)f\F] dzdr<K{5t)^ 



and 



sup sup 



E 



(i~s)t<^<t J (i-S)t Jo 
0<y<l 



E' 



Dr,zu{iy,y) 



dzdr < K{Sty- 



(3.15) 



(3.16) 



P-almost surely. 



Proof: We will only deal with the proof of (I3.15p . since (I3.16p may be checked using exactly 
the same arguments. 

Let us first invoke the linear equation (13. lift satisfied by the Malliavin derivative Du{u,v), 

for {y, v) E [(1 — 6)t, t] x (0, 1), and then take the square L^-norm on [(1 — 5)t, t] x [0, 1]: 

! [ \Dr,zu{iy,v)\^ dzdr <2a'^ f [ \G^-riv, z)\^ dzdr 

J{l-5)t Jo ' J{l-5)t Jo 

+ 2/ [ (f f G,.s{v,y)b\u{s,y))Dr,Ms,y)dyds] dzdr. (3.17) 

J{l~5)tJ0 \Jr Jo / 

By (13.10ft . the first term in the right-hand side of (13.17ft can be bounded by ^J^a'^{6t)^. For 

the second one, we apply Holder's inequality, the fact that b' is bounded and Fubini's theorem, 
so that we end up with 



\Dr^z'iJ'{T^,v)\'^ dzdr < 



il-S)t Jo 



(l-5)t Jo 



\Gu-s{v,y)\' 



{l-5)t Jo 



\Dr,z'u{s,y)\'^ dzdr \ dyds. 
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Taking the conditional expectation E[ ■ \F] and using again (I3.10p we obtain 

sup [ [ E[\Dr,MP^v)\^\^] dzdr < if^a^dt)^ 

-6)t<p<v J(l^S)t Jo ' V TT 



0<'«<1 



+ 



/II 2 



6t 



sup 

(l-5)t \{i-6)t<T<s J {l-S)t Jo 
0<y<l 



E [\Dr,M'r,y)\'^\F] dzdr 



ds. 



If we set 



sup [ [ E[\DrMP,v)\^\F] dzdr, u e[il-6)t,t\, 
{i-s)t<p<^ Jn~5)t Jo 

0<v<l 



then we have seen that 

^5,t(i^) <K{6t)^ + KSt 



1 

^5,t(s)^^= ds, z/ G [(1 - 5)t,t], a.s. 

(l~S)t - s 



Now we can conclude by applying Gronwall's lemma [5^, Lemma 15]. □ 

Proof of Proposition \3.2\: We first recall that the Malliavin derivative of u(i',v), {iy,v) G 
[0,r] X [0,1], satisfies that Dr^zu{i',v) > 0, for all {r, z) G [0,T] x [0,1], a.s. This is because 
the Malliavin derivative solves the linear parabolic equation (13. lift . Let us deal with the proof 
of (13.141) in two steps: 

Step 1: Ttie lower bound. Fix 6 G (0, 1] and let us first derive the lower bound in (I3.14p . 
Since the Malliavin derivative of u{t,x) is non-negative, formula (13. 13^ yields 



9{F) > [ e-'E E'( f [ DrMt, ^) (aS?^)) dzdr 

Jo L \J(l-5)t Jo ^ ^ 



F 



de. 



By Equation (13.111) . we can decompose the right-hand side of the above inequality in a sum 
of four terms: 



Ao{t, X] 5) = 



{i-&)t Jo 



\Gt-r{x, z)\'^ dzdr, 



Ai{t,x;5) 



-a 



Jo 

X 



Gt-r{x, Z 
ft /.I 



dzdr. 



A2{t, x; 5) =a 



Gt-ri^i z) 



X E 



E' 



(l-5)t Jo 
t rl 




Gt-s{,x,y)h\u{s,y)) (Dr,zu{s,y)) dyds 



r JO 



F 



dzdrdO, 
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A3{t,x;6) 



E 



X E' 



(l-5)t JO 
t rl 



t /.I 




Gt-s{,x, y)b'{u{s, y))Dr,zu{s, y) dyds 



r JO 




Gt-s{x,y)b'{u{s,y)) {Dr,zu{s,y)) dyds 



F 



dzdrdO. 



'r JQ 

First we notice that, by (I3.10p : 

Ao{t,x;S) >a^c,idt)'^, (3.18) 

Thus we can write 

g{F) >a^c^iSt)^^ - \A,{t, x; 6) + A2it, x; S) + As{t, x; S)l (3.19) 
so that we will need to obtain upper bounds for the terms \Ai(t, x;6)\, i = 1,2, 3. 

We apply Fubini's theorem for the conditional expectation, the boundedness of b' , Cauchy- 
Schwarz's inequality and the bound (I3.10p . so we have the following estimates: 

\Ai{t,x;6)\ 



< C 



t 

{l-S)t Jo 

t 

{1-S)t 
t 



\Gt-rix, z)\'^ dzdr 

mGt-s{x,y)E[\Dr,Ms,u)\ \F] dyds 
\Gt-s{x,y)\^E [\Dr,Ms,y)\^ \F] dyds ) dzdr 



dzdr 



1 / rt /•! 



(1-5)1 Jo \Jr Jo 




(1-5)1 Jo 



\Gt-s{x,y)\' 



{l-5)t Jo 



E [|Dj,,z'u(s, y)p |F] dzdr j dyds 



<C{5t)\ sup [ ! E[\Dr^zu{s,y)\'^\F]dzdr 

\(i-S)t<8<t J {i^S)t Jo 



0<i/<l 



At this point we are in position to apply (13.151) in Lemma [231 Therefore 

\Ai{t,x]5)\ < G {5t)i, a.s., 



(3.20) 



where G 



^a^b'lU 



0F(27r)l 



In order to get a bound for \ A2{t, x;6)\, one can use analogous arguments as for \ Ai{t, x; 5) \ 
but applying (13.161) instead of (I3.15p in Lemma (13.30 . Hence, one obtains 



\A2{t,x]6)\ < G (5t)3, a.s. 



(3.21) 
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Let us finally estimate \A3{t, x; S)\. For this, we apply Fubini's theorem, the fact that b' 
is bounded, Cauchy-Schwarz inequality with respect to dzdrdP\pdP' and we finally invoke 
Lemma 13.31 



\A^{t,x-5)\ 



< C 



1 rt 




{l-5)tJ0 J{l-5)tJ0 



Gt-s{x,y)Gt--s{x,y) 



E 



< G 



{1-S)t Jo 

oo 



E' (\Dr,zu{s,y)Dr,Ms, : 



dzdr j dydsdydsdO 



1 rt 




{l~5)tJ0 J{l~5)tJ0 



(1-5)1 Jo 
t 



E 



Gt-s{x,y)Gt^s{x,y) 
1 

V 

dzdr 



{l-5)t Jo 



E 



E' {\Dr,Ms,y)\' 



< G{5t)\ 



Gt-s{x,y)dyds 



dydsdydsdO 



{i-S)t Jo 



where G 



1 ^2||L/||2 



(3.22) 

The very last estimate in (13.221) has been obtained after applying 

Cauchy-Schwarz's inequality and the bound (I3.10p . Eventually, plugging the bounds (I3.20p - 
(!T^ in (ITT^ we have 



where Ci is a positive constant depending on a and ||6'||oo- Hence, if we assume that 5 <l/\j 
then we can write 



T' 



g{F) > (T^c^{5t)2 - 2ci(5t)3 > a'cj^ - 2ci5-iT^ . 



-5 „3 \ 1 



-5_3 . 



It only remains to observe that the quantity a'^c^Si — 2ci6iTi is strictly positive whenever 6 
is sufficiently small, namely 6 G (0,(5o), with 

(5n = 1 A — A 



T T V 2ci 

Thus, the lower bound in (I3.14p has been then proved. 

Step 2: The upper bound. The upper estimation in (I3.14p is almost an immediate conse- 
quence of the computations which we have just performed for the lower bound. More precisely, 
according to (13.130 and the considerations in the first part of the proof, we have the following: 



g{F)<J2\Mt,x;l)\, 



(3.23) 



i=0 
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where we notice that we have substituted 5 by 1 in Ai{t, x; S), i = 0, 1, 2, 3. We have already 
seen that \Ai{t, x;l)\ < Ct* , for i = 1, 2, and \A^{t, x]5)\ < Ct^ , so we just need to bound 
|y4o(t, x; 1)|, which follows directly from ( I3.10p . Thus 



for a constant C > 0. Therefore g{F^ < C2t'^ , where the constant C2 > depends on a, 
and T. □ 

We are now in position to prove the main result of this section: 

Proof of Theorem \3.2i The random variable F = u{t,x) — E{u{t, x)) is centered, belongs 
to 2 and, by Proposition [O it holds that < CiU < g{F), for all t G (0,T]. We apply 
then [19j Theorem 3.1] and we obtain that the probability density p : M ^ M of F is given by 

E\uit,x)-E{uit,x))\ fry. 
for almost every 2 G M. Thus, the density p of the random variable u{t, x) satisfies 

2g(z-E(u{t,x))) ^\ Jo g{y) ) 

In order to conclude the proof, we only need to use the bounds (13.140 in the above expression 



4 The stochastic heat equation in 

in this section we are interested in the stochastic heat equation in with an additive Gaussian 
noise which is white in time and it has a spatially homogeneous correlation in space. We aim 
to find out sufficient conditions on the drift term and the noise's spatial correlation ensuring 
that the density of the solution admits Gaussian-type estimations. The fact that we deal with 
a SPDE in with a non-trivial spatial correlation makes the analysis in this case much more 
involved in comparison to the one-dimensional setting in Section [31 

4.1 SPDEs with spatially homogeneous noise 

Let us consider the following stochastic parabolic Cauchy problem in W^: 

— {t,x)-Au{t,x) = b{u{t,x))+aW{t,x), (t, x) G [0, T] x M"^, (4.25) 

tJL 

where T > 0, A denotes the Laplacian operator in R'^, b : M ^ M is a function with 
bounded derivative, and suppose that we are given an initial condition of the form 

u{0,x) = uo{x), x G R'^, 
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with : M*^ ^ R measurable and bounded. 

The random perturbation W (formally) stands for a Gaussian noise which is white in 
time and with some spatially homogeneous correlation. More precisely, on a complete prob- 
ability space P), we consider a family of mean zero Gaussian random variables W = 
{W{(p), (f G C^(R"'+i)}, where C^(M°'+i) denotes the space of infinitely differentiable func- 
tions with compact support, with covariance functional 

E{Wiip)Wi^j))= / iipit)*^sit))ix)Aidx)dt, ip,^eC^iR''^'), (4.26) 
Jo Jr'^ 

where ^psit,x) := tlj{t,—x) and A is a non-negative and non-negative definite tempered 
measure. By |25l Chapter VII, Theoreme XVII], A is symmetric and there exists a non- 
negative tempered measure /i whose Fourier transform is A. That is, by definition of the 
Fourier transform on the space iS'(]R'^) of tempered distributions, for all belonging to the 
space i5(]R'^) of rapidly decreasing C°° functions. 



and there is an integer m > 1 such that 

/ (1 + |eP)-XrfO < oo. 
Jm.''- 

The measure /i is usually called the spectral measure of the noise W . In particular, we have 







Notice that we have used the symbol "*" for the standard convolution in 



Example 4.1 Usual examples of spatial correlations are given by A{dx) = f{x)dx, where f 
is a non-negative continuous function on M.'^ \ {0} which is integrable in a neighborhood ofO; 
for instance, one can take f to be a Riesz kernel fe{x) = for < e < d. The space-time 

white noise would correspond to consider f equal to the Dirac delta on the origin. In this 
latter case, the spectral measure /i is the Lebesgue measure on W^. 

We denote by 7i the completion of the Schwartz space i5(M"') endowed with the semi-inner 
product 



)h= (0i*(02).)(x)A(rfx)= / 



01,02 ^ and associated semi-norm || ■ H^^. Notice that 7i is a Hilbert space that may 

contain distributions (see [5, Example 6]). Set Tir '■= L^([0, T]; 7-^). 

Then, it turns out that the Gaussian noise W can be naturally extended to TYt, so that 
we obtain a family {W{h)^ h G H-t} of centered Gaussian random variables such that 

E{W{h^)W{h2)) = {h^,h2)n, = r{h{t),h2{t))ndt, huheHr. 

Jo 
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The details of that extension can be found, for instance, in p. 805, or [23], p. 3. Then, 
if we set 

Wt{g) := te[Q,T],geH, 

we obtain a cylindrical Wiener process {Wt{g), t G [Q,T],g E H} on the Hilbert space H. 
That is, for any t E [0,T] and g eH, Wt{g) is a mean zero Gaussian random variable and 

EiWtig,)Ws{g2)) = (t A s){g,, g^U, s,tE[0,T], gi,g2en. 

We denote by the a-algebra generated by {Ws{g), s E [0,t],g E H} and the P-null sets. 
As it has been explained in PH Section 3], one can construct (real-valued) stochastic integrals 
of predictable processes in L^(fi x [0, T]; H) with respect to the cylindrical Wiener process W; 
notice that we are making an abuse of notation since W denoted the Gaussian family given 
at the beginning. The resulting stochastic integral turns out to extend Walsh's integration 
theory [26J and it is equivalent, in some particular situations, to Dalang's stochastic integral 
set up in [5]. The stochastic integrals appearing throughout this section will be understood 
as integrals with respect to the cylindrical Wiener process W. 

We are now in position to rigorously define the mild solution of Equation (14.250 : a square 
integrable {jF^j-adapted random field {u{t,x), (t, x) E [0,T] x R'^} solves Equation (14.251) 

if, for any (t, x) E [0,T] x M"', 



Gt{x - y)uo{y) dy + [ [ G. 

Jo Jwi 



u{t,x)= Gt{x - y)uo{y) dy + / Gt-s{x - y)b{u{s,y) dyds 



+ a / Gt-s{x-y)W{ds,dy). (4.27) 
Jo Jr"^ 

Here, the function Gt{x), t > 0, x G M°', denotes the fundamental solution associated to the 
heat equation in W^, that is the centered Gaussian kernel of variance 2t: 

1 _^ 
GAx) = e « . 

General existence and uniqueness results for Equation (14.270 may be found in |5]. More 
precisely, it turns out that sufficient conditions for existence and uniqueness of solution to 
Equation (14.270 are the following: h is Lipschitz, uq is measurable and bounded, and the 
noise's spatial correlation is related with the fundamental solution G through the condition 

/ / \J'Gt{i)\'' ii{d^)dt<+oo. (4.28) 
Jo 

As it has been shown in [5", Example 8], condition (14.280 holds if and only if 

^-^^^{dO < +00, (4.29) 
and this integrability condition will be assumed to be satisfied in the remainder of this section. 
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Let us turn now to the question whether the probability law of the solution at any point 
has a density. The Gaussian setting in which we apply the Malliavin calculus machinery is 
determined by the Gaussian family {W{h), h G Ht} given before. 

Using this framework, it is a consequence of |2I1 Theorem 5.2] that, if the drift coefficient 
b E has a bounded and Lipschitz continuous derivative and condition (14.290 is fulfilled. 



then the solution u(t,x) to Equation (14.271) . at any point {t,x) E {0,T] x W'-, is differentiable 
in the Malliavin sense, that is u(t,x) E D^'^, and its law has a density with respect to the 
Lebesgue measure. At this point, we should mention that in [21j all the results are proved 
in the case where the noise's correlation is given by a function /, that is A(dx) = f{x)dx. 
However, the extension of those results to a general tempered measure A is straightforward. 
See also the works [T5], [211 ^O'' I'slated results with slightly stronger conditions on the 
spectral measure. 

It will be of much importance for us the equation satisfied be the Malliavin derivative of 
u{t,x). Indeed, see either [15j, [23J or [21j, the Malliavin derivative Du{t,x) takes values in 
the Hilbert space Ht and satisfies the following linear parabolic equation: 

Du{t,x)=aGt-.{x-^)+ [ [ Gt-s{x-y)b'{u{s,y))Du{s,y)dyds, (4.30) 

Jo JR-i 

where 'V stands for the 7Y-variable. Moreover, one proves that 

sup sup E{\\Du{t,x)\\l^J < +00. (4.31) 

te[o,r] xeRd 

Equation (14.301) may be interpreted in the following sense: for any r E [0, t), Dru(t, x) satisfies 
the equation in H 

Dru{t,x) = aGt^rix - -k) + / / Gt-^six - y)b'{u{s,y))Dru{s,y) dyds. (4.32) 

Jo Jr'' 

The integral on the right-hand side of (I4.32p is understood as a TY-valued pathwise integral. 
Before going on with our analysis, let us briefly describe how this integral is rigorously defined. 
Let {cj, j > 1} be a complete orthonormal system of H. Then, using the properties of G, 
the boundedness of b' and (14.311) . the 7i-valued integral 

Zt= / Gt-s{x -y)b'{u{s,y))Dru{s,y) dyds 
Jo Jr'^ 

can be defined through its components 

aGt~s{x - y)b'{u{s,y)){Dru{s,y),ej)'H dyds, j >1 

with respect to the basis {ej, j > 1}, and those latter integrals takes values in M. Moreover, 
one can obtain an upper bound for the square moment of (for the general setting see, for 
instance, [22], p. 24): 

E{\X,\^)<C f [ \G,^s{x -y)\'\b'{u{s,y))\'E{\\DMs,y)rn) dyds. (4.33) 

Jo JR-i 
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Let us go back now to Equation (14.321) . Since Gt-r{x — ★) : R'^ — > M defines a function 
(indeed in 5(M'^)), this implies that either Dru{t,x) and the integral in the right-hand side 
of (I4.32p define elements of H which are functions in z. Therefore, for any fixed {t,x) E 
(0,T] X W^, we can state that the Malliavin derivative satisfies, for all {r, z) G [0,t) x Mf^: 

Dr,zu{t,x) = aGt-r{x - z) + / Gt-s{x -y)b'{u{s,y))Dr,^u{s,y) dyds. (4.34) 

Jo JR'* 

A crucial consequence of this fact is that, as in the case of the one-dimensional stochastic 
heat equation with boundary conditions, the Malliavin derivative Dr^zu{t,x) is non-negative, 

for all (r,z) G [0,t) X M'^, a.s. 

We will need a slightly stronger condition on the spectral measure /i than (14.291) . Namely, 
consider the following hypothesis: 



Hypothesis H^: There exists r] G (0, 1) such that 

1 



/i(dO < +00- (4-35) 



Then, one can prove the following estimates (see [T5l Lemma 3.1]): 
Lemma 4.2 Assume that the spectral measure fi satisfies J^^ j^^Kd^) < +oo. 

1. LetT > 0. Then, there exists a constant ki > such that, for all t e [0,T], 

kit< [ f \J^G siOl^ fi{dOds. 
Jo Jr-^ 

The constant ki depends on T and, indeed, it converges to zero as T tends to infinity. 

2. Suppose that Hypothesis Hrf holds. Then, there exists a constant ^2 > such that, for 
all t > 0, 

ft 



/ / \J^Gs{0\^ ^l{dOds<k2t^ (4.36) 
Jo Jm.<^ 



for all (3 e (0,1-7]]. 



Remark 4.3 It is worth mentioning that the integrability condition ^4.28\) was sufficient for 
us to prove the existence of density for the solution u(t, x), at any point {t, x) G (0, T] x W'- 
(see j22 Theorem 5.2]). However, as it will be made clearer in Section \4~^ we will really 
need upper bounds of the form ( 14.36]) in order to obtain Gaussian estimates for the density of 

u{t, x). 
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4.2 Gaussian estimates for the density of the solution 

Let us consider T > and {u(t, x), (t, x) G [0, T] x R'^} the unique mild solution to Equation 
( I4.27p . This section is devoted to proof the following result: 

Theorem 4.4 Fix t E (0,T] and x E M'^. Suppose that Hypothesis is satisfied for some 
7] E (0, |) and that the coefficient b is of class and has a bounded Lipschitz continuous 
derivative. Then, the random variable u(t, x) has a density p with respect to Lebesgue measure 
which satisfies the following: for almost every 2; G M, 

, , E\u(t,x)—m\ f (z — mY^ 

where m = E{u{t,x)) and Ci, C2 sre positive constants depending on a, ||6'||oo- V snd T. 

Theorem 14.41 will be a consequence of [121 Theorem 3.1] and the following proposition. As 
we have done in Section 13.21 we use the notation F = u{t,x) — E{u{t,x)) and we remind 
that we will need to find almost sure lower and upper bounds for the random variable g{F), 
where 

poo 

g{F) = / e-'E 
Jo 

Proposition 4.5 Fix T > 0. Assume that Hypothesis H^ holds for some rj E (0, |) and the 
coefficient h is of class and has a bounded Lipschitz continuous derivative. There exist 
positive constants Ci,C2 such that, 

Cit < g{F) < Cst^-", a.s., (4.37) 

for alltE (0,r]. 

In order to prove Proposition 14.51 we will need the following technical lemma, which plays 
the role of Lemma ESI in our standing setting. 

Lemma 4.6 Let t > and assume that Hypothesis Hrj holds. Then, there exists a positive 
constant C depending on a, ||b'||oo 3nd the constant k2 in ( I4.36|) , such that, for all 5 E (0, 1], 

sup E 

{1-S)t<v<t 

and 

sup sup E 

0>1 {1-S)t<u<t 

for any /5 G (0, 1 - ?]]. 



E \u{t, x) 



exp 



m] 



Cit 



E' ({Du{t,x),Du{t,x)) 



Hi 



de. 



\DrU(u, 



n dr 



Ja-5)t 



F 



< C{6tf, a.s. 



(4.38) 



E' 



\DrU{u, 



n-5)t 



n dr 



F 



< C{Sty, a.s. 



(4.39) 
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Proof: It is very similar to that of Lemma 13.31 In fact, owing to Equation (14.320 . we have, 

for any {u,v) G [(1 - 6)t,t] x M"^, 



{l-5)t 



\Dru{u,v)\\li dr < 2(7^ 



\\G^^r{v dr 



+nh'\\Lm 



{l-S)t 



(i-&)t 



D,.u{s,y)\\l^dr] dyds, (4.40) 



where we have applied Minkowski's and Cauchy-Schwarz's inequalities. By (14.361) in Lemma 



{l-5)t 



- *)||^ dr = r / 1^(01' KdOdr < H^tY 
Jo Jr"^ 



for all /? G (0, 1 — ?7]. Therefore, plugging this bound in (14. 40 ft and taking conditional expec- 
tation, we obtain: 



E 



iDrUiu, v) IlL dr 



{1-S)t 



F 



+ 



1 



/l|2 



2V27r 



( 



(l-5)t 



sup E 

(l-«)t<T<S 



v 



' y)\VH dr 



\ 



1 



Vi^ - s 



ds, a.s. 



As we have done in the proof of Lemma 13.31 we are now in position to apply Gronwall's lemma 
IS Lemma 15]. Hence (14.38ft is proved. The estimation (14.39ft can be checked using exactly 
the same arguments. □ 

Proof of Proposition 14.51 ' The framework of the proof is similar to that of Proposition 13.21 
in Section [T21 However, computations here will be slightly more involved since we are working 
in a Hilbert-space-valued setting determined by Ht = L'^{[0,T];'H) . Let us first deal with the 
lower bound in (I4.37p : 

Step 1: The lower bound. Recall that F = u{t,x) — E(u{t,x)) and the random variable 
g{F) can be written as 



9{F) 



e-'E 



e-'E 



E' ({Du{t,x),Du{t,x))nT) \F 



dO 



E' \ I {Dru{t, x), Dru{t, x))n dr 



F 



d9, 



where Du(t,x) denotes the shifted random variable {Du{t, x)){e ^uo + \/l — e~'^^io') 
According to Equation (14.30ft . for any 6 G (0, 1] we have the decomposition 

g{F) > a^Bo{t,x] 6) - \Bi{t,x] 6) + B2{t, x; 6) + Bs^t, x; 6) \ , 



(4.41) 



where 
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Bo{t,x;S)= / \\Gt^r{x - dr, 

l{l~5)t 



Bi{t,x; 6) 
= aE 



t I rt 

I Gt-r{x - -k 
{l-S)t \ J{l-5)t 



, / Gt^s{x - y)b'{u{s,y))Dru{s,y) dyds) dr 

J(i-s)t Jr'' I -w 



) dr 


F 


n 





B2{t,x;S) =a / e'^E 
Jo 



{l-5)t 



E'\ I {Gt-r{x-i.) 

(1-5)* 





F 







dd, 



Bsit,x;S)=l e~'E 
'o 



E'{ I (I I Gt-s{x-y)b'{u{s,y))Dru{s,y)dyds 

{1-S)t \J{l-5)t 



{l-5)t 



) dr^ 


F 


de 









By part 1 in Lemma |4~^ notice first that 

Boit,x;5)= [ [ \J^Gs{0? ii{dC)ds>ki5t. 
Jo Jr^ 



(4.42) 



Concerning the second term Bi{t,x;6), we can apply Cauchy-Schwarz and Mikowski's 
inequalities, so that we obtain 



t \ 2 

2 



\B^{t,x;5)\<G {6t)-2 [ / - ★)||^dr 



X E 



{1-S)t 



{1-S)t 



Gt-s{x - y)b'{u{s, y))Dru{s, y) dyds 



2 






dr 


F 


) 


H 







-5t 



<Ci6t)-^ / / \J^GriO\ f^idOdr 
\Jo Jr'^ 

X I r / \Gt-s{x-y)\'E 

'{l-5)tJRd lJ{l-S)t 



|D^M(s,?/)||^rfr 



dyds 



Thus, by (I4.36p . Lemma |42] and the fact that 

\Gt^s{x-y)\Uyds<CiSt)k^ 



(l-5)t 



(4.43) 
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we have: 

\Bi{t,x;6)\ <C{6tf+-^, (4.44) 

for all /3 e {0,1 — r]]. The term \B2{t, x; 5)\ can be treated in the same way as we have just 
done for \Bi{t, x; 6)\. Namely, one proves that 

\B2{t,x;6)\ < C {5tY^ ( f \\Gt^r{x - ^)\\l,dr^ ' 

\J{l~5)t 



X / e-' 
Jo 

< C 

X e 
Jo 



E 



E' 



[1-S)t 



{i-&)t 





" \ 




V 


s{x - y)b'{u{s, y))Dru{s, y) dyds 


dr j 


F 


dO 




n ) 








\Gt-six-y)\'E 



(l-5)t 



E' { I \\Druis, 

{l-5)t 



dr 



F 



Taking into account (14.361 ). Lemma 1431 and (I4.43p . we also get 



\B2{t,x-5)\ < C{5t) 



P+7 



dyds \ dO. 



(4.45) 



for all (3 e {0,1-7]]). 

Eventually, in order to deal with the term iB^^t, x; 5)\, we mainly apply Cauchy-Schwarz's 



inequality with respect to E 



E' In- 



F 



so that we will take advantage of the 



computations which we have performed so far. More precisely, we have 

B3it,x;6) 

ft ft f 

Gt-s{x - y)h'{u{s, y))Dru{s, y) dyds 



< E 



POO 

X e 
Jo 



(\-&)t 
E 



2 






dr 


F 


) 


H 







E' 



{1-S)t 



{i-&)t 





" \ 






s{x - y)b'{u{s, y))Dru{s, y) dyds 


dr j 


F 


de 




H ) 







The two terms in the right-hand side of the above inequality already appeared in the analysis 
of Bi{t,x;5) and B2{t,x;S), respectively, and each of them may be bounded, up to some 
constant, by ((5t)2 + 4. Therefore 



\Bsit,x; 6)\<Gi6tf^K 
for any /3 e (0, 1 - r]]. Estimations ( ITTO . (l4l^ and (^M)-<^M) yield 

g{F) > (r^h5t - ci (^(5t)^+t + , 

where ci depends on a, ||&'||oo and Hence, \f 6 < 1 A ^ and /9 G (^, 1 — rj], 

9{F) > t (^a^k.d - 2ci(5^+tt/3-i j > t (^a^k,6 - 2ci6^+-^T^~-^^ . 



(4.46) 
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Observe that the quantity 

defines a positive constant whenever 6 E (O.Sq), with 

^ 1 1 f(j^ki\^- 
5q = 1 A - A - ' * 



T T V 2ci 

Therefore, we obtain the desired lower bound in (14.370 



Step 2: The upper bound. The upper in (14.370 is an almost immediate consequence of 
the computations in the Step 1 and the decomposition 

3 

g{F)<Y,mt,x-l)\. (4.47) 

1=0 

Indeed, we have already found upper bounds for the last three terms on the right-hand side of 
( I4.47p . On the other hand, observe that (I4.36p yields 



for all /3 e (0, 1 - 77]. This bound, together with (I4.44ft - (l4.46p in the case 6 = 1, implies 



9{F) <c{tl' + t^+i + t^+i) < C2t^ 



for all 13 E (0, 1 — 77], where the constant C2 depends on T. Therefore we conclude the proof. 

□ 

The proof of Theorem 14.41 can be finished as in the case of Theorem 13.11 □ 

5 Gaussian bounds for the stochastic wave equation at 
small time 

The main objective here is to extend the results in Section |4]to a stochastic wave equation in 
space dimension c? < 3 and controlled by the spatially homogeneous Gaussian noise considered 
there. The intrinsic properties of the differential operator driving the equation will not allow 
us to obtain optimal results for all time parameter T > 0, even if we assume that the noise's 
space correlation satisfies Hypothesis H,^, an assumption which is slightly stronger than the 
one needed to have existence of density for the corresponding mild solution (see [21j). 

5.1 The stochastic wave equation in dimension <i = 1, 2, 3 

We consider here the same setting as in Section |4] but for the stochastic wave equation in 
spatial dimension d < 3: 

— {t,x)-/\u{t,x) = h{u{t,x)) + aW{t,x), (t, x) G [0, T] x M'^, (5.48) 

LJL 
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where T > 0, 6 : R — M is a function with bounded derivative, and suppose that we are 
given initial conditions of the form 

u{0, x) = Uo{x), — (0, x) = fo(x), X G M'^, 

with uo,vo : M*^ ^ M measurable and bounded functions such that uq is of class C^(]R'^) 
and has a bounded derivative Vmq- The random perturbation W corresponds to the spatially 
homogeneous Gaussian noise described in the previous Section |4~T1 We recall that fi denotes 
the corresponding spectral measure and {J-'t, t > 0} the filtration defined by the cylindrical 
Wiener process associated to the noise W. 

The mild solution of Equation (I5.48p is given by a {jF^j-adapted process {u{t,x), (t, x) G 
[0,T] X R'^} such that, for all {t,x) G (0,r] x R'^, 

u{t, x)= [ Vo{x- y)Ti{dy) + ^( [ Uo{x - y)Tf{dy)] 

jRd VJRd / 

+ b{u{s,x-y)TU{dy)ds + a f [ TU{x - y)W{ds, dy), (5.49) 

Jo Jwi Jo JRd 

where Tf, t > 0, denotes the fundamental solution of the wave equation in dimension d = 



1,2,3: 

(a;) = ^l{|x|<t}, 

The element at stands for the surface measure on the three-dimensional sphere of radius t. In 
particular, for each t, Tf has compact support and, in the case d = 3, is no more a function 
but measure on R^. It is important to remark that only in these cases, T^, d = 1,2,3, defines 
a non-negative measure. Existence and uniqueness of mild solution to Equation (15.490 is a 
consequence of the results in whenever the space correlation satisfies 



1 + 



-^(dO < +00. (5.50) 



We also point out that the stochastic integral in the right-hand side of ( 15.490 is a well- 
defined integral of a deterministic element in Ht with respect to the cylindrical Wiener process 
associated to the noise (see Lemma 3.2 and Example 4.2 in |21J). 

For all dimensions c? > 1, we have a unified expression for the Fourier transform of F^: 

, _ sin(27rt|g|) 

Using this fact and assuming that ( 15.501) holds, one proves the following lemma (see [13J, 
Lemmas 5.4.1 and 5.4.3): 



21 



Lemma 5.1 For anyt > it holds that 

with some positive constants ci, C2 > 0. 

Thus, if we assume that J^^ i+]fp < (15.511 ) yields 

di{tAt^)< [ [ \Tr'',{0\^^i{dOds<d2{t + t^), 

Jo jRd 

for all t >0, with some positive constants di,d2. In particular, for t E [0, 1) we have 

dit' < [ f |J^rf(0|V(t^0^^s<4t. (5.52) 
Jo Jr^ 

However, under Hypothesis H^ (see (I4.35P ) one can get a slightly sharper upper estimation 
(see [23, Lemma 3]): 

Lemma 5.2 Let T > and assume that Hypothesis holds. Then 

f I \TTi{i)\^ ii{di)ds<d,t^-^\ (5.53) 
Jo Jr'^ 

for all te [o,r]. 

Eventually, if = 1, 2, 3, explicit computations yield that, for any t > 0, 

/ / Vi{dy)ds < Ct^ , (5.54) 

Jo Jr'^ 

where C is a positive constant that only depends on d. 

If b belongs to and it has a Lipschitz continuous bounded derivative, then the solution 
u{t,x), at any {t,x) G (0,r] x R'^, belongs to D^'^ and its Malliavin derivative, as a random 
variable taking values in Ht = L'^{[0,T];H), satisfies 

Du{t,x) = aV^^Xx -i^) + [ [ b'{u{s,x -y))Du{s,x-y)rt^^{dy)ds, (5.55) 

Jo iiRd 

where 'V stands for the 7Y-variable (see [21, Proposition 5.1]). This linear equation is 
understood in x [0,T];7Y) and let us remark that crr^_.(x — ★) is a well-defined element 
in Ht (see PH Lemma 3.2]). Moreover, under the standing hypothesis, the random variable 
u{t,x), for {t,x) e (0,T] X R"', has a density with respect to the Lebesgue measure (see [2T1 
Theorem 5.2]). Of course, the Gaussian setting here is the same as the one that has been 
considered in Section IH 



22 



5.2 Gaussian estimates of the density at small time 

For T > 0, consider the unique mild solution {u{t,x), {t,x) E [0,T] x M'^} to Equation 
( I5.48p . In this section, we will prove that the density p : R'' — > R of u{t,x) has lower and 
upper Gaussian bounds whenever T is small, where this essentially means that T < 1. The 
main result is the following: 

Theorem 5.3 Suppose that Hypothesis /-/^ is satisfied and that the coefficient b is of class 
and has a bounded Lipschitz continuous derivative. Then, there exists Tq E (0, 1] such that 
the following statement is satisfied: for any T E (0, Tq) and (t, x) E (0, T] x R"*, the random 
variable u{t, x) has a density p with respect to Lebesgue measure such that, for almost every 

zeR, 

E\u(t,x) — m\ { (2 — m)^ 1 , , E\u{t,x) — m\ { (z — mY 



where m = E{u{t, x)) and Ci, C2 are positive constants depending on a, ||&'||oo. Tq and rj. 

Remark 5.4 In the case of the stochastic heat equation presented in Section\^ we have been 
able to obtain Gaussian upper and lower bounds for any T > 0, while here we restrict our 
analysis to small T. As we will precisely point out in the next Proposition \5.5[ that difference 
is due to the fact that the Malliavin derivative of the solution to the stochastic wave equation 
does not need to be a non-negative function. 

The statement of Theorem 15.31 is an immediate consequence of [121 Theorem 3.1] and the 
following proposition. For t > and x E R"^, set F = u{t, x) — E{u{t, x)) and we remind that 
we will need to find almost sure lower and upper bounds for the random variable g{F), where 

g{F) = J^ e-^E E' (^{Du{t,x),TMt^x))nT^\F de. (5.56) 

Proposition 5.5 Assume that Hypothesis Hr^ holds. There exist Tq E (0, 1] and positive 
constants Ci, C2 such that, for any T E (0, Tq), 

Cit^ < g{F) < C2t^~^\ a.s. (5.57) 

for all tE [0,T]. 

Proof: It follows the same lines as the proof of Proposition |431 so that we will only point out 
the main steps. More precisely, we observe first that the Malliavin derivative Du(t,x) solves 
the linear equation ( 15.551) with a non-negative initial condition but driven by a hyperbolic 
operator. Thus, in comparison with the stochastic heat equation, Dr^zu{t,x) does not need to 
be non-negative as a function of (r, z); indeed, in the case d = 3, even the Malliavin derivative 
does not need to be a function. Hence, in order to deal with the lower bound of g{F) (see 
(I5.56P ). we will not be able to restrict the integral with respect to dr on a small time interval 
as we have done in the proof of Proposition 14.51 This is the reason why we will be forced to 
consider T < 1. 
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In fact, by (15.551) . we are only able to consider the decomposition 

giF) > Doit) - {\D,it)\ + \D2it)\ + \D,m : 



(5.58) 



where 



Di{t) = aE 



D2{t) 



Doit) = a' f WrUix - ^)r^ dr, 
Jo 



Ttrix-^), 




Jr'^ 



\ dr 


F 


n 





E' { I /T^^ 



TUix-^) 



^ j^b\u{7^y)){Dru{^-y)^ Tt^{dy)ds'^ dr^ 



de. 



e~'E 



E' i I ( [ [ b'{u{s, X - y))Dru{s, x - y) Vt^{dy)ds 

\-'0 J9.<i 



F 



de. 



nb'{u{s,x - y)) (Dru{s,x - y)) T'}_^{dy)ds\ dr 
'■•^ ^ ^ I n ^ 

By the lower bound in (15.521) . we have 

Doit) > dit^. (5.59) 
Concerning the term Di, we can argue as follows: 

\Dr{t)\<c( f\\TU{x-^)\\l,dr 



X E 



< Ct-^ E 



5-2?7 

< ct-^ 



< ct^-'^\ 




b\u{s, X - y))Dru{s, x - y) T'l^Xdy)ds 
f(f I \\D,u{s,x-y)\\nTU{dy)ds\ dr 

Jo \Jo JR<* / 



Jm.'^ 
t / ft 



2 






dr 


F 


) 


n 








E 



JR'i 



\Dru{s,x- y)\\lidr 



IJO 



TUidy)ds 



(5.60) 



where we have used (15.531) . (I5.54p and (15.631) in Lemma ESI below. 
Using similar arguments one proves that 

\D2it)\<Ct'-^^. 



(5.61) 



24 



The analysis of |L'3(T)| can also be performed by following the calculations above to obtain 
(I5.60p . so that we end up with 

\D-i{t)\<Cf~^''. (5.62) 
Plugging the estimates (ESS]) -((5221) in (ESHD yields 

9{F) > d.t' - C3 (t^-^" + t^-^") , 

for all t E [0,T], where C3 is a positive constant depending on a, ||fo'||oo and rj. Hence, if 
T < 1 we have 

9iF) > f {d, - 2c,T'-'^) , 
and the quantity Ci := di — 2oiT'^^'^^ is strictly positive whenever T < Tq, where 



Tn = 1 A 



dl \ 2-2" 

2^ 



Therefore, we have proved the lower bound in (I5.57p . 

The upper bound in (I5.57p is an immediate consequence of what we have done so far and 
(I5.53p . because 

3 

9{F)<j2mt)\<c2t'-''^- 



i=Q 



□ 



In the proof of Proposition 15. 5[ we have applied the following technical lemma, whose 
proof is very similar to that of Lemma 14.61 

Lemma 5.6 Let t > and assume that Hypothesis Hn holds. Then, there exists a positive 
constant K depending on a, \\h'\\oo snd the constant d^ in Lemma lS^ such that 



sup E 

0<s<t 



\DMs,y)\\'^ dr 



(5.63) 



and 



sup sup E 

e>i o<s<t 



E' 



\DMs,y)\\y^ dr 



F 



(5.64) 
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